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Is dynamo in tachocline or in convection zone 1 ? 
That is the question! 
Whether 'tis nobler in the mind to fit 
Parameters of model to reproduce the observations 
Or to take the theory and direct simulations 
And try to answer basic questions: 
Why magnetic field is generated? 

And why it drifts equatorward in course of solar cycle? 
'tis a consummation devotely to be wished! 

take-off Shakespeare 



Abstract. We give a short introduction to the subject and review advances in understanding the 
basic ingredients of the mean-field dynamo theory. The discussion includes the recent analytic 
and numerical work in developments for the mean electromotive force of the turbulent flows and 
magnetic field, the nonlinear effects of the magnetic helicity, the non-local generation effects 
in the dynamo. We give an example of the mean-field solar dynamo model that incorporates 
the fairly complete expressions for the mean-electromotive force, the subsurface shear layer and 
the conservation of the total helicity. The model is used to shed light on the issues in the solar 
dynamo and on the future development of this field of research. 
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1. Introduction 



The mean-field magnetohydrodynamic presents one of the most powerful tools for 
exploring the nature of the large-scale magnetic activity in cosmic bodies flMoffatt||l"978 



Parker|1979}|Krause fc Radler|1980l ). It is widely believed that magnetic field generation 



there is governed by interplay between turbulent motions of electrically conductive fluids 
and global rotation. The dynamo theory studies the evolution of the magnetic field which 
is govern by the induction equation: 

dB 

Ik 



V x (U x B — 77 V x B) 



where B is the magnetic field induction vector, U is the velocity field of the plasma and 
rj is the molecular magnetic diffusivity. 



2. The mean electromotive force 

The aim of this section is to briefly outline the basic equations and methods of 
the mean-field magnetohydrodynamic (MHD). The general framework of the mean-field 
MHD can be introduced as follows. In the turbulent media, it is feasible to decompose 
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the fields on the mean and fluctuated parts, e.g., B = B + b, U = U + u. Hereafter, 
everywhere, we use the small letters for the fluctuating part of the fields and capital 
letters with a bar above for the mean fields. Let's define the typical spatial, L,£ , and 
temporal, T, r c , variation scales for the mean and the fluctuated parts of the fields. It is 
a typical situation for the astrophysical system when the flows and magnetic fields are 

UTq 

strongly turbulent, i.e., the Strouhal number St = — — ~ 1 and the Reynolds number 

ul 

R e = — S> 1, where v is the molecular viscosity. Assuming the validity of the Reynolds 



rules (Monin & Yaglom 19751 and averaging the induction equation over the ensemble 



of the fluctuating fields we get the mean-field dynamo equation: 

— = V x (UxB + S) , £=u^Tb (2.1) 

There are two contributions here. In a perfectly conducting fluid, when the magnetic 

ul 

Reynolds number R m = — 3> 1, the magnetic flux is frozen into fluid. Then, the first 

V 

term can be interpreted as the defection, stretching and compression (or expansion) of 
the magnetic field by mean flow, because of V x (UxB) = — (U ■ V) B + (B • V) U — 
B (V • U) . The effect of the turbulence is represented by the mean electromotive force 
£—u x b. It is possible to analyze the general structure of the £ using the assumption 
about the scale separation in the turbulence L,T 3> £, t c and the transformation sym- 



metry properties of the basic physical quantities (Radler 1969 Krause & Radler 1980 
Brandenburg et al.||2012[ ): 



S = (a + j) o B - fj o (V x B) + k o (VB) + o ( - ) , (2.2) 

where (VB)j. .j — | (^iBj + VjSj), the kinetic coefficients a, 7, fj, k are tensors and 
the symbol o marks the tensor product. 

The kinetic coefficients may depend on the global factors, which determine the large- 
scale properties of the astrophysical system, for example, the global rotation angular 
velocity i~2, the large-scale vorticity W = VxU, the stratification parameters like Vp, 
Vu 2 , and the global constraints, like, magnetic helicity conservation. For th e simplest 



case when a ~ ao^j, i) ~ r\T§ij and 7 ~ — Sij n V n and k is neglected we get(|Krause & 



Radler||l980"l ): 



Et = aoBi + (V (p) xB)- t (VxB), (2.3) 

— (p) 

where uq is the magnitude of the a effect, V is the turbulent pumping velocity and 
t\t is the isotropic turbulent diffusivity. 

To calculate the kinetic coefficients we use the equations which govern the evolution 
of the fluctuating magnetic and velocity fields. For example, taking into account the 
effects of the global rotation and shear for the incompressible turbulent flows, we get the 



equations as follows (see, e.g., 



Radler et aL]|2003 



3 

Brandenburg fe Subramanian||2005 ) 



9b 

dt 



= Vx(uxB + Uxb)+ 7]V 2 b + 0, 



diij 
dt 



+ 2(n x u). = -Vi [p 



(2.4) 
(2.5) 



/<+& 



where 25, J stand for the nonlinear contributions of fluctuating fields, p is the fluctuating 
pressure, f2 is the angular velocity responsible for the Coriolis force, f is the random 
force driving the turbulence. 

Using Eqs fl2.4|2.5D , S can be calculated analytically by different ways. The first-order 
smoothing, also known as the second order correlation approximation (SOCA) uses the 
condition min (Rm, St) <C 1 and neglect the nonlinear contributions in Eqs( 2.4|2~5 ) ( |Mof-| 



fatt 1978 Krause & Radler 1980). The r-approximations was introduced to take into 
account the nonlinear effects of the second order correlations. It is claimed to be valid 
for Rm, Re 3> 1 and for the developed turbulence in equilibrium state. In this case we 
solve the equations for the second order correlations and replace the third-order correla- 
tions of the fluctuating parameters, with the second order relaxation terms (see details 
in |Kleeorin et al.| [19961 |Blackman fc Field|[2002l |Radler et al] [20031 |BrandenburglZ 
Subramanian||2005 ) . This approximation is based on questionable assumptions ( Radler 



& Rheinhardt 2007), e.g., it is assumed that the second-order correlations do not vary 
significantly on the time scale of r c . This assumption is consistent with scale separation 
between the mean and fluctuating quantities in the mean-field magnetohydrodynamic. 
The reader can find a comprehensive discussion of the r -approximation in the above 
cited papers. 



The path integral approach use the ideas from the stochastic calculus (Dittrich et al. 
1984| IZel'dovich et al.||1984[ ). This approach is vali d for the case Re > 1, Rm < 1. The 
reader can find the relevant examples in papers by Kleeorin & Rogachevskii ( 1999 ) and 



by Rogachevskii et al. (20111. The every analytic al method calculation of £ has to use 
the assumptions about the background turbulence which would exist in the absence of 
the large-scale magnetic fields and flows (e.g., global rotation and shear). For the nu- 
merical solution it is equivalent to definition of the stochastic driving force f . Despite 
the Eqs( 2.4|2~5 ) is widely applied to the dynamo in the Sun and the late- type stars, 
these equations describe the forced isothermal turbulence rather than turbulent convec- 



tion. The latter also can be treated analytically using the r-approximation (Kleeorin & 
Rogachevskii]|2003[ ). 



The mean-electromotive force can be estimated by the direct numerical solution (DNS) 



of the equations like Eqs( 


2.4|2.5 


) (e.g., 


Brandenburg|2001 


Kapyla & Brandenburg||2007| 


Brandenburg et al.||2008a 


Livermore et al.||2010 


Tobias et al.||2011 


) or the similar ones 



for the turbulent convection by the so-called "impose- field" method (e.g., Ossendrijver 
et al.|2001[|2002[) or the so-called "test-field" method (|Schrinner et al.|2005||Kapyla et al. 



2008[|2009|}Rheinhardt fc Brandenburg|2010|[Schrinner|2011||Schrinner et al.|2012[ ). The 
global simulations of the geo- and stellar dynamos also can be used to extract the mean- 



field dynamo coefficients from simulations (see the above cited papers and Racine et al. 
2011||Brown et al.||201l| . 
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3. Mean-field phenomena in the solar magnetic activity 

There is a wide range of the magnetic activity phenomena which are observed on the 
Sun and the others astrophysical systems. Here, I restrict myself with consideration to 
the solar magnetic activity. The observation of the solar magnetic fields shows that for 
the spatial scales the basic assumption behind the Eqs( 2.1|2~2 | is not fulfilled (see, e.g., 
the review by Sami Solanki in this volume). The given theory is not able to capture 
self-consistently and simultaneously the origin of the large-scale sunspots butterfly dia- 
grams and the emergence of the separate sunspots. Though the both phenomena can be 
analyzed separately using the mean-field MHD framework, see, e.g., the application of 
the theory to the sunspot decay problem) Riidige r fc Kitchatinov|2 000). It is known that 
the large-scale temporal-spatial (e.g., time-latitude) patterns, such as the sunspots but- 
terfly diagrams, can be detected for the much smaller scales phenomena, like ephemere 
regions ( pakarov fc Makarova||1996| |Harvey||2000| |Makarov et al.||2004[ ). Therefore, even 
the scale-separation assumption is not valid for the solar conditions we can consider the 
large-scale organization of the magnetic activity phenomena as a manifestation of the 
large-scale magnetic fields generated somewhere in the deep convection zone. This idea 
is commonly adopted in the mean-field dynamo theory. 

The dynamo theory isolates of the details of processes, which are responsible for the 
emergence of the magnetic activity features at the surface, and studies the evolution of the 
large-scale magnetic field govern by the dynamo equations Eqs( 2.1|2T2 ). It is suggested 
that the toroidal part of the large-scale field forms sunspots and organize the magnetic 
phenomena inside the Sun. The large-scale poloidal field goes out of the Sun and governs 
the solar corona. Thus, the key questions for the theory are the origin of the large-scale 
magnetic activity spatial-temporal patterns, the phase relation between activity of the 
poloidal and toroidal components, what defines the solar cycle period and magnitude 
etc. Another portion of the problems which could be studied using the same framework 
is related to the statistical properties of the large-scale spatial organization of the small- 
scale magnetic fields and motions in the solar convection zone with the operating dynamo 
(see, e.g, Parnell et al. 2009 Stenflo & Kosovichev 2012 and review by Jan Stenflo in 



this proceedings). 



The basic idea for the solar dynamo action was developed by Parker ( 1955 ) . He sug- 



gested that the poloidal field of the Sun is stretched to the toroidal component by the 
differential rotation (f2 effect) and the cyclonic motions (a effect) return the part of the 
toroidal magnetic field energy back to the poloidal component. This is the so-called afi 
scenario. While the turbulent diffusion is not presented in the title, it is equally important 
( Parker|[l979| . The mean electromotive force given by Eq(|2.3|) fits to this scenario. The 
fi effect can be well understood because the helioseismology provide the data about the 
distribution of the rotation inside the convection zone and beneath. The big uncertainty 
is about how the poloidal field of the Sun is generated. There is an ongoing debate about 



a number of problems connected with the a effect and afl dynamos (see, e.g., Rudiger 
fc Hollerbach|[2004 Brandenburg fc Subramanian|2005 ). For instance, the period of the 
solar activity cycle poses a problem. Namely, for mixing-length estimates of the turbulent 
magnetic diffusivity in the convection zone and dynamo action distributed over the whole 
convection zone, the obtained cycle periods are generally much shorter than the observed 
22 yr period of the activity cycle. For thin-layer dynamos, the situation becomes even 
worse. 

The expression for the mean electromotive force £ contains a number of dynamo effects 
that may complement the a effect or may be an alternative to it. These effects are due 
to a large-scale current, global rotation and (or) the large-scale shear flow. The given 
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dynamo effects are usualfy associated with the 17 x J effect (Radler|l969l and the shear- 
current effect W x J where, W = V x U ( Rogachevskii fc Kleeorin|2003[). I n fact, these 
effects contribute to the antisymmetric parts of fj in the £ given by Eq(2.2|. The reader 
can find the explicit expressions for them in ( Radler et al.|20 03 Rogachevskii & Kleeorin 



2003] |Pipin|[2008| ) . [Pipm" fc Seehafer| ( [2009] ) and |Pipin fc kosovichev| ( |2011a| ) found that 



the inclusion of the additional turbulent induction effects increases the period of the 
dynamo and brings the large-scale toroidal field closer to the equator, thus improving 
the agreement of the models with the observations. Also, in the models the large-scale 
current dynamo effect produces less overlapping cycles than dynamo models with a effect 
alone. The symmetric part of k (see, Eq. 
diffusivity (see, e.g., 



Kichatinov et al. 1994 



2.2) contribute to the anisotropic turbulent 
Rg,dler et aL|2003| |Pipin||2008[ ) It is rather 



important for the dynamo wave propagation inside the convection zone (Kitchatinov 
20CJ2] |Kichatlnov1[2053| 



The DNS dynamo experiments support the existence of the dynamo effects induced 
by the large-scale current and global rotation ( Schrinner et al.|2005 Kapylii et al. 2008). 
It was found that we have to account the complete expression of the £ ( see, Eq. 2.2 ) to 
reproduce the simulations of the global dynamo action ( Schrinner et al.||2005 Schrinner 



2011 Morin et al.||2011 ) and evolution of the large-scale fields in the convective rotating 



turbulent flows ( |Kapyla et al.|2009||Brandenburg et al.|2012[ ». The aim of the DNS is to 

simulate the dynamo action in the cosmic bodies and we are still on the way to reproduce 
the basic properties of the large-scale dynamo for the Sun. It was shown that the mean- 
field MHD framework is useful for the analysis of the results obtained in simulations (see, 
also, Racine et al.||2011). 



4. The magnetic helicity issue 

The properties of the symmetry transformation of the £ suggest (Krause & Radler 



19801 that the a effect is pseudoscalar (lacks the mirror symmetr y) which is related to 

T~c 

the kinetic helicity of the small-scale flows, i.e., ao — u • Vxu. 

3 



Pouquet et al. 



|1975| 

showed that the a effect is produced not only by kinetic helicity but also by the current 
t c ( — b • V X b 

3 ( U ' VXU 2^T" 

as resistance of magnetic fields against to twist by helical motions. It leads to the concept 
of the catastrophic quenching of the a effect by the generated large-scale magne tic field 

a (0) 



helicity, and it is ao = 



. The latter effect can be interpreted 



It was found that ao (B) = 



1 + R m (B/B eq ) 



(Kleeorin & Rogachevskii 



of R m S> 1, the a effect is quickly saturated for the large-scale magnetic 



1999). In case 
strength 



field 



that is much below the equipartition value B eq ~ y p/iqu 2 . The result was confirmed by 
the DNS(Ossendrijver et al. 2001). The catastrophic quenching (CQ) is related to the 
dynamical quenching of the a effect. It is based on conservation of the magnetic helicity, 
X = a-b (a is fluctuating part of the vector potential) and the relation between the 
current and magnetic helicities he = b • V X b ~ x/^ 2 , which is valid for the isotropic 
turbulence! Moffatt 1978). The evolution equation for \ can be obtained from equations 



that governs a and b, it reads as follows (Kleeorin & Rogachevskii 1999 Subramanian 
|fc Brandenburg||2004| : 



dx 
dt 



-2 £-B) 



V • T x - 77B • J, 



(4.1) 



6 



V.V. Pipin 1 



where we introduce the helicity fluxes T x = a x u x B ax (uxb). The helicity fluxes 
are capable to alleviate the catastrophic quenching (CQ). The first example was given 
for the galactic dynamo model (Kleeorin et al. 2000). The calculations and the DNS 



shows the existence of the of the several kind of the helicity fluxes. Part of them have the 
turbulent origin, e.g., the fluxes due to the anisotropy of the turbulent flows (Kleeorin & 
Rogachevskii 1999), the fluxes due to the large-scale shear (Vishniac & Cho 2001) and 
the diffusive fluxes ( |Mitra et al.|2010" ). Another kind of the helicity fluxes, which are not 
mentioned in the Eq.(4.1 ), are related to the large-scale flows, e.g., meridional circulation 
and outflows due to the solar wind ( Mitra et al.|2011 ). Generally, it was found that the 
diffusive fluxes, which are ~ %Vx, where rj x is the turbulent diffusivity of the magnetic 
helicity, work robustly in the mean-field dynamo models but it requires r\ x > tjt to reach 
\B\ > 0.lB eq . 

Another possibility to alleviate the catastrophic quenching is related with the non-local 



formulation of the mean-electromotive force ( 


Brandenburg & Sokoloff 2002 


Brandenburg 


et al. 


2008b 


)■ 


Kitchatinov & Olemskoy 


(20111 found that the nonlocal a effect and the 



burg & Kapyla, (2007) show that the result can depend on the model design. Nonlocal 



formulation of the mean-field MHD concept suggests a possibility to solve the problem re- 



lated to the dynamo period in the mean-field models. Rheinhardt & Brandenburg (2012 1 
showed that corrections for nonlocal effect in the mean-electromotive force can be taken 
into account with the partial equation like (l + r c dt + £ 2 A) = £„• , where£„-°' ) is the 



local version of the mean-electromotive force given by Eq.(2.2l. 

We have to notice, the solar dynamo is an open system, where the large-scale mag- 
netic fields escape from the dynamo region to the outer atmosphere. For the vacuum 
boundary conditions, which are widely used in the solar dynamo models, the magnetic 
field escapes freely from the solar convection zone, and nothing prevents the magnetic 
helicity accompanying the large-scale magnetic field to escape the dynamo region. Thus, 
the magnetic helicity conservation should not pose an issue for the solar type dynamos. 
Recently, Hubbard & Brandenburg (2012) revisited the CQ concept and showed that 



for the shearing dynamos the Eq.(4.1) produces the nonphysical fluxes of the magnetic 



helicity over the spatial scales. Hubbard & Brandenburg (2012) suggested to cure the 



situation starting from the global conservation law for the magnetic helicity, 



j t J {a-b + A-B}d^ = -77 J {B- J + b-j}dV- J 



V-T x dV 



(4.2) 



where integration is done over the volume that comprises the ensemble of the small-scale 
fields. We assume that J- x is the diffusive flux of the total helicity which is resulted from 
the turbulent motions. Ignoring the effect of the meridional circulation we write the local 
version of the Eq.fl42j) as follows ( |Hubbard fc Brandenburg||2012[ ): 



d t ab = —dt (A • B) 



X 



Rr, 



? ? B J V-T x 



Note, that the large-scale helicity is govern by: 

dt (A • B) = 25 • B + V • ((£ x A) - A 



(U x B)) 



(4.3) 



(4.4) 



Therefore, Eq s.p~l)) and p~3| ) differ by the second part of Eq. p~4| . [Hubbard fe Bran 
denburg (2012) found that the V • (£ X A) cures the problem the nonphysical fluxes of 



the magnetic helicity in shearing systems. Another term V • (A x (U x B)) contains the 
transport of the large-scale magnetic helicity by the large-scale flow. It was found that 
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the dynamos with the dynamical quenching govern by the Eq.(4.3) does not suffer from 
the catastrophic quenching issue. 



5. The dynamo shaped by the subsurface shear layer 

Most of the solar dynamo models suggest that the toroidal magnetic field that emerges 
on the surface and forms sunspots is generated near the bottom of the convection zone, 



in the tachoclinc or just beneath it in a convection overshoot layer, (see, e.g., Ruediger 
fc Brandenburg|[T995l |Bonanno et~aLl[2002{ |Rempel||2006l |Guerrero fc de GouveiaTM 



Pino 



2008 



2009). However, an attention was drawn to a number of theoretical and 



observational problems concerning the deep-seated dynamo models ( Brandenburg||2005 



2006 1 . We proposed (see Pipin fc Kosovichev|2011c|b I a solar dynamo model distributed 



in the bulk of the convection zone with toroidal magnetic-field flux concentrated in a 
near-surface layer. In this section I would like to illustrate the theoretical profiles for the 
principal parts of the mean-electromotive force for the linear and nonlinear case in the 
mean-field large-scale dynamo and compare them with the results of the DNS. 



We study the mean-field dynamo equation Eq.(2.2) for the axisymmetric magnetic 
field B — e^B + V X (Ae^/ (r sin 9)) , where is the polar angle. The expression for 
the mean electromotive force £ is given by Pipin (2008): 



Si = 



Bj — i]ijkS7jBk- 



(5.1) 



The tensor onj represents the a-effect. It includes hydrodynamic and magnetic helicity 
contributions, 



C a sin 2 9a { " ) 



a 



(AT) 



= s t 



{zvt (/ 10 a) ( 

i {** (ft' 

3 VT {(e.A^ 



+ e^e 



e • A 



(p) 



(5.2) 
(5.3) 



(a) 
6 



e,:A 



(u) 



ej A 



/s (a) }, 



where the hydrodynamic part of the a-effect is defined by a\, , = Vlogp quanti- 
fies the density stratification, A^ = C„Vlog \ Vt^ quantifies the turbulent diffusivity 
variation, and e = O / | O | is a unit vector along the axis of rotation. The turbulent 
pumping, 7^ , depends on mean density and turbulent diffusivity stratification, and on 
the Coriolis number ft* = 2t c VLq where r c is the typical convective turnover time and 51 
is the global angular velocity. We introduce the parameter C v for the Aw to take into 



account the results from the DNS (Ossendrijver et al. 2001 Ka/pyla et al. 20091 which 



show that the a effect is saturated to the 
the profiles below. Following the results of 



bottom 



Pipin 



of th e convection zone. We will show 



(|2008|), 7^ is expressed as follows: 



7£ ) =3^{/J" ) A«+/« (e-AW)e n } 



3Wi 



ien A&>, (5.4) 



3^ 



The effect of turbulent diffusivity, which is anisotropic due to the Coriolis force, is given 
by: 



fid) 
J2 



Hjk 



2/i 



(a) 



njk 



(<0 



(5.5) 
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Figure 1. Left, the profiles of the a effects components for the 9 — 45°. Right, the profiles of 
the background turbulent diffusivity C v tit, the isotropic, rf 1 ^ , a nd anisotropic, rf A \ pa rts of 
the magnetic diffusivity and Q. x J effect (also known as S effect (Radlcr 1969; S tix|1976 1. 



The functions fj^f^j in Eqs(5.2 5.4|5.5) depend on the Coriolis number. They can be 
found in |Pipin| ffiiUty (see also, |Pipin fe Kosoviche'v1|2011b| |Pipin fc Sokoloff||2011[ ). In 

b 2 

the model, the parameter e = =, which measures the ratio between magnetic and 

_ Popu 2 

kinetic energies of the fluctuations in the background turbulence, is assumed to be equal 



to 1. In our models we use the solar convection zone model computed by Stix (2002). 

The mixing-length is defined as I = Q!mlt|A( p )| , where A( p ) = Vlogp quantifies 
the pressure variation, and «mlt = 2. The turbulent diffusivity is parameterized in the 



form, 



- r r7 (0) 



where rffi 



is the characteristic mixing-length turbulent 



3 /w 

diffusivity, I is the typical correlation length of the turbulence, C v is a constant to 
control the efficiency of large-scale magnetic field dragging by the turbulent flow. Also, we 
modify the mixing-length turbulent diffusivity by factor f s (r) = 1 + exp (100 (r ov — r)), 
r ov — O.72i?0 to get the saturation of the turbulent parameters to the bottom of the 
convection zone. The latter is suggested by the DNS. The results do not change very much 
if we apply AW = C„Vlog [Wt^ with C v ^ 0.5. For the greater C v we get the steady 
non-oscillating dynamo concentrated to the bottom of the convection zone. I would like 
to stress that the purpose to introduce the additional parameters like C v — 0.5 and 
f s (r) is to get the distribution of the a effect closer to the result obtained in the DNS. 
The bottom of the integration domain is rt, = 0.71-Rq and the top of the integration 
domain is r e = O.99i?0. The choice of parameters in the dynamo is justified by our 
previous studies ( Pipin fc Seehafer|2009 Pipin fc Kosovichev|2011a ), where it was shown 
that solar- types dynamos can be obtained for C Q /C W > 2. In those papers we find the 
approximate threshold to be C a ss 0.03 for a given diffusivity dilution factor of C v — 0.1. 
Figure [T] shows the radial profiles for the principal components of the mean electromotive 
force, which are essential for our model. They are in the qualitative agreement with the 



results of the DNS obtained by Ossendrijver et al. (2001) and Kapyla et al. (2009). The 



contribution of small-scale magnetic helicity \ = a 'b to the a-effect is defined as 



AM) _ Aa) 



XTc 



Hop£ 2 ^pl 2 



(5.6) 



The nonlinear feedback of the large-scale magnetic field to the a-effect is described by 
a dynamical quenching due to the constraint of magnetic helicity conservation given by 



Eq.(4.3 ). For the illustration purpose we use the realistic value for the magnetic Reynolds 



Figure 2. Snapshots of the magnetic field and helicity evolution inside the North segment of 
the convection zone. Top panel panel shows the field lines of the poloidal component of the mean 
magnetic field and the toroidal magnetic field (varies ±700G) by gray scale density plot. The 
bottom panel shows the large-scale (density plot) and small-scale magnetic helicity (contours) 
distributions. 



number R m — 10 6 . We matched the potential field outside and the perfect conductivity at 
the bottom boundary with the standard boundary conditions. For the magnetic helicity 



the number of the possibilities can be used (Guerrero et al. 2010 Mitra et al. 2010). 
We employ V r (x + A • Bj = at the top of the domain and % = at the bottom 
of the convection zone to show that for the a-quenching formalism which is based on 
the Eq.(4.3), the boundary conditions determine the dynamics of the dynamo wave at 
the near surface layer. To evolve the Eq.(4.3) we have to define the large-scale vector 
potential for each time-step. For the axisymmetric large-scale magnetic fields where the 
vector-potential is 



A = e A P + vT = 



r sm ( 



-A + re r T. 



(5.7) 



The toroidal part of the vector potential is governed by the dynamo equations. The 
poloidal part of the vector potential can be restored from equation V x (rT) = e^B. 

Figure [2] shows the snapshots of the magnetic field and magnetic helicity (large- and 
small-scale) evolution in the North segment of the solar convection zone. The Figure 
shows the drift of the dynamo waves which are related to the large-scale toroidal and 
poloidal fields towards the equator and towards the pole, respectively. The distributions 
of the large- and small-scale magnetic helicities show one to one correspondence in sign. 
This is in agreement with Eq.(4.3). It is seen that the negative sign of the magnetic 
helicity follows to the dynamo wave of the toroidal magnetic field. This can be related 
to the so-called "current helicity hcmispcric rule" which is suggested by the observations 



(ISeehafer 


1990 Zhang et al. 


2010). The origin of the helicity rule has been extensively 


studied in the dynamo theory 


(Fisher et al.|1999 


Choudhuri et al.|2004 Kuzanyan et al. 


2006; Sokoloff et al.||2006[ |Pevtsov & Longcope||2007 


Pipin & Kosovichev||2011b 


Zhang 



et al. 2012). 



The time- latitude diagrams for the dynamo model are shown in Figure [3] The results 
arc in qualitative agreement with observations. We show the dynamical a effect as well. 
The model shows that with the given boundary conditions the a effect increases and has 
positive maximum at the growing phase of the cycle and it decreases, having the negative 
minimum at the decaying phase of the cycle. The variations of the radial profiles for the 
a effect and the small-scale-magnetic helicity are shown in Figure [4] The saddle in the 
the a effect profile is resulted from the given boundary conditions and distribution of 
the large-scale magnetic field near the top of the solar convection zone. The latter is 
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Figure 3. Top, the time - latitude variations of the toroidal field near the surface, r — 0.95Rq, 
(contours ±300G) and the radial magnetic field at the surface (density plot). Bottom, the same 
for the the toroidal field and the a effect (a^ component) (density plot). 




0.05 



0.75 0.80 0.85 0.30 0.95 1.00 50 100 150 201) fa JKI 350 100 150 

DISTANCE, [r/R a ] TIME, [YR] 

Figure 4. Left panel shows variations of the a effect (a^ component) profiles with the cycle at 
colatitude — 45° , the kinetic part of the a effect is shown by the solid line. Right panel shows 
the ratio of the maximum of the toroidal field strength and the equipartition value with time. . 



determined by the dynamo boundary conditions and by the subsurface shear layer. It 
was found that the saddle disappears in case of V r x = at the top. The highly non- 
uniform radial distribution of the a effect (at least for the near-equatorial latitudes) was 
found in the global LES (large-scale eddy) simulation of the dynamo action in the recent 



paper by Racine et al. (20111 



The results obtained with the given dynamo model can be summarized as follows. 
In the model the dynamo wave propagates through the convection zone and shaped by 
the subsurface shear layer. The model incorporates the fairly complete expressions for 
the mean-electromotive force, including the anisotropic turbulent pumping and magnetic 
diffusivity, the turbulent generation of the poloidal and the toroidal fields due to a 2 and 
Q x J effects and the dynamical a quenching due to magnetic helicity.We demonstrate 
that the conservation of the total helicity given in the form of Eq.(4.3l alleviates the 
catastrophic quenching (see Figure |4| for the solar type dynamos that operates on the 
base of the afl mechanism. The further development can be related to the nonlocal 
framework for the mean-electromotive force. 
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6. Conclusions 

Summarizing the main topics of this review, which I also consider as a principal ad- 
vances in the mean-held dynamo theory, I conclude as follows. Firstly, the knowledge of 
the mean-held coefficients is decisive in order to analyze and model dynamo action in a 
variety astrophysical bodies. It was found that agreement between the direct numerical 
simulation and mean-held models is good only if a large number of mean-held coefficients 



are involved, which contribute to a, 7, 77, k (see, Eq 2.2 1. One obvious reason for this is 
that the spatial scale-separation is not valid for many astrophysical cases. Secondly, the 
issues related to the catastrophic quenching of the dynamo can be alleviated by different 
ways. However, the approach, which is based on conservation of the total helicity, see 



Eq.(4.3l, suggests just a more than the way around the catastrophic quenching issue. It 
was proved that the magnetic helicity conservation does not pose an issue for the solar 
type dynamos. The new formalism can be used to study the origin of the different kind 
of anomalies in solar magnetic activity, which are observed and related to evolution of 
the large-scale toroidal held, the surface distribution of the magnetic helicity at different 
phase of the solar cycle etc. These points have not been discussed in this review and 
remain an open held for the future work. 
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